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We investigate on exotic tetraquark hadrons of the kind [cs] [cs] by computing their spectrum and 
decay modes within a constituent diquark-antidiquark model. We also compare these predictions 
with the present experimental knowledge. 



Introduction. In the last few years we are witnessing the discovery of a number of new narrow hadronic 
resonances with charm which do not match standard quark-antiquark interpretation, thereby named exotic 
"q^" ■ hadrons. This has called for alternative interpretations of their inner structure. One of the possible explanations 
is that these particles are loosely bound molecules of open charm mesons Another possibility is that new 
aggregation patterns of quarks in matter are possible. We follow the suggestion by Jaffe and Wilczek of having 
C^l ' di-quarks as building blocks 2]. Light diquarks have been object of several lattice studies. Recently we have 
!— i . studied the hypothesis of treating the diquark as a fundamental field in large N Chromodynamics in two 
Oh' dimensions Q finding some qualitative features of the spectra which are common to real ones. The idea that 
, the colored diquark can be handled as a constituent building block is at the core of the approach taken in this 
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Motivated by our former study on the interpretation of the V(4260) resonance [4|, we analyze the possibility 
r— !■ of a spectroscopy of particles with hidden strangeness and charm embodied in diquark-antidiquark structures 
; of the kind [cs] [cs] , where q = [cs] is a 3 C diquark. We predict the mass spectrum for these states and discuss 
i*- which might be their prominent decay modes on the basis of quark rearrangements in the qq system. The 
O i' mass spectrum is computed using the non-relativistic spin-spin interactions Hamiltonian, supposed to remove 
the degeneracy among the various [cs] [cs] states with assigned spins and anguar momenta 
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The couplings are inversely proportional to quark masses: we will incorporate this dependency in the color- 
magnetic moments k. The Hamiltonian {!]) describes contact interactions. For this reason we expect that 
allowing a relative orbital angular momentum between the diquarks will decrease or switch-off the spin-spin 
interactions between quarks and antiquarks. We will take into account this effect in the determination of the 
CS| ' spectrum of the negative parity particles. 

The model. Adopting the approach discussed at length in [j| and [(J we will determine the mass spectrum 
of [cs][cs] hadrons by diagonalization of the following non-relativistic effective Hamiltonian for a [9132] [Q1Q2] 
diquark-antidiquark hadron 
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H = 2m q + H^f + H^f + H SL + H L (2) 



where 



-^SS^ = 2k ?i9 2 ' + ' ^2 ) + 2Hq 1 q 1 S qi ■ Sq i + ^q 2 q 2 Sq 2 ■ Sq 2 (3) 

represent the chromomagnetic interactions between quarks in the tetraquark system whereas the spin-orbit and 
orbital contributions are given by 

H SL = 2A{S^ -L + S^-L) 
L(L+1) 

H L = B 1 2 ' (4) 

respectively. The symbol 5 q represents the total spin of the diquark q = [91(72]- A,B are coefficients to be 
determined by data. To diagonalize H we need to estimate the diquark mass m q and the coupling constants 
k; then we have to specify the states having assigned J PC quantum numbers and find their masses. The 
[cs] diquark mass is simply estimated by ?7ir cs i = m\q S ] + m c — m q = 1955 MeV, where m c q are constituent 
masses and mw is obtained by comparison with light scalar mesons data 0- We have used mwi = 590 MeV, 
m c = 1670 MeV, m q = 305 MeV. 
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As shown in |5j, quark-antiquark spin-spin couplings are estimated by comparison with the mass spectra of 
ordinary qq mesons. A calculation made along these lines provides us with the chromomagnetic couplings for 
qq color singlets {n qq ) 1 . We obtain n cs = 25 MeV, (k cs -)i = 72 MeV, (k 8I )i = 121 MeV, (k c5 )i = 59 MeV, 
A = 22 MeV, B = 495 MeV. 

On the other hand, the couplings n qq in Eq. ([3]) are not necessarily in the singlet channel as those estimated 
since octet couplings (ks) are also possible between quarks and antiquarks in a qq system [T3|. The octet 
couplings are estimated with the aid of the one-gluon exchange model as follows. For the diquark attraction 
is in the 3-color channel, we can write q l = [cs]' 1 := e 1; ' k CjS k , neglecting spin. i,j,k are color indices in the 
fundamental representation of SU(3). Then the color neutral hadron is 

[cs}[cs] = e ijk e ij <k>(cjS k )(c j 's s ') = { Cj c>){s k s k ) - ( Cj s j )(s k c k ) (5) 

We then use the following SU (N) identity for the Lie algebra generators 
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J] A&Ag, = 2 5 u S jk - -6 l3 6 kl (6) 



where N is the number of colors. A color octet (N=3) qq state can be written as q^Xf.q^ , and consequently 



^X^o>)(s k X a kl s l ) = J2 A^cVW = 2 

a 

This allows to extract from §5§ the octet term as follows: 

[cs][cs] = {c j c l ){s k s k ) 



(c l s t )(s k c k ) - ^(c lC l )(s k s k ) 



(7) 



±(?X? j cj)(s k Xl l s l ) + ±(c j ci)(s k s k ) 



2 -{cfc>){s k - S k ) - ^(c i Xi J c j )(s k Xl l s l ) (8) 



This formula gives information about the relative weights of a singlet and an octet color state in a diquark- 
antidiquark tetraquark. We have three colors running in the sum qc' whereas = 1, ...,8 in cfA^-c 3 . Therefore 
the probability of finding (projecting onto) a particular qq pair in color singlet, for example to find cc in the color 
singlet state Cj(? , is half the probability of finding the same pair in color octet cX a c as 3 x 2/3 = 1/2(8 x 1/2). 
We write then: 
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K C c{[ cs ][ cs }) = g( K cc)l + g(Kcc)8 



(9) 



where (« c g)i have been reported above. For the determination of the quantity (k c5 )s we have to resort to the one- 
gluon exchange model. In this model we assume that the coupling K qq iq)(R.), R being the color representation 
of the two quark system can be considered proportional to the product of the two color charges in the vertices 
of a qq{q) — > qq{q) diagram in which one gluon is exchanged between quarks. No distance dependency is taken 
into account. 

Then, if k cE (R) is the weight of the quark-antiquark interaction, writing the above mentioned product of 
color charges in terms of the SU(3) Casimir operators we have 



Kce (R) ~ (C( 2 >(R) - C^(3) - C^(3)) 

We recall that C (2) (R) = 0, 4/3, 4/3, 3 as R = 1, 3, 3, 8. Then it is immediately found that 
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(Kcc)i ~ —Z (K c s) f 
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(«cc)j 



(10) 
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Finally, from Eq. ©, we have 



Kcc([cs][cs]) = -(«cc)l 



(12) 



Now that we know its input parameters, we are ready to diagonalize the Hamiltonian ^j. We label the 
particle states by using the notation 1 5 q , 5q ; , J) where S^q is the total spin of the diquark-antidiquark 
system. States are organized in order to have definite J PC quantum numbers. For negative parity ones a unit 
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of relative angular momentum between the diquark and the antidiquark is required (£ qq — 1). Altogether we 
have 

a-Two positive parity states states with J PC = ++ 

|0 ++ ) 1 = |0 cs ,0^;0 qq! J = 0) |0++) 2 = |l cs ,l es ;O qq , J = 0) (13) 

b-Three states with J = and negative parity (L qq = 1 required) 

\A) = |l cs ,0 Es -;l qq , J = 0) \B) = |O cs ,l s ;l qq , J = 0) \C) = |l cs , l a ; l qq , J = 0) (14) 

State |C) is even under charge conjugation. Taking symmetric and antisymmetric combinations of states \A) 
and \B) we obtain a C-odd and a C-even state respectively; therefore we have two states with J PC = h 

\0- + ) 1 = ^(\A)-\B)) \0- + ) 2 = \C) (15) 



and one state with J PC = 



\0-) = -j=(\A) + \B)) (16) 



c-Three states with J — I and positive parity 

| J D) = |l CS) gs -;l qq ,J=l) |£) = |0 cs ,le 5 ;l qq ,J=l) |F) = |l cs ,l sg ;l qq ,J=l) (17) 

\F) is an eigenvector under charge conjugation, with negative eigenvalue. Operating on \D) and \E) in the same 
way as for states \A) and \B) we obtain the J PC = 1 ++ state 

|l ++ > = -^ (|£>> + |i?» (18) 



and the J PC = 1+ ones 



\1 + -) 1 = ^(\D)-\E)) \l + -) 2 = \F) (19) 



d-Six states with J = 1 and negative parity. To start with consider the following two 

\G) = |l cs ,0 es -;l qq , J= 1) \H) = \0 cs ,l- cS ;l m ,J =1) (20) 

differing from \D) and \E) as we have L qq = 1 here. When symmetrized and antisymmetrized these give the 
combinations 

- ^= (|G> - W)) II — >! = ^= (IG) + |if)) (21) 

Moreover we have the following four charge conjugation eigenstates 

|l-+) 2 = |l cs ,l Es -;l qq , J= 1) |1— ) 2 = |0 CS) a -;0 qq ,J=l) 

|1 — > 3 -|l cs ,l es ;0 qq ,J=l) |1— ) 4 = |l cs ,l^;2 qq ,J=l) (22) 

The action of the spin operators in Eq. ([3]) on the states here listed is independent of the specific L qq value. 
Let us write 

|<5 q , 5 q ; i5 qq , J) = |c Ts,c Ts]S^,J) (23) 

where the T can be T° = 1/V2 <j 2 and T l = 1/V2 a 2 <j % for spin and spin 1, respectively. The numerical factors 
are chosen in such a way to preserve the normalization Tr[(r Q )^T' 5 ] = 5 a/3 . Then the action of a spin-spin 
interaction operator, e.g. S c ■ S s , on the generic state in Eq. (|23p is: 

(S c ■ S s )\c T Ts, c T rs; S qq , J) = ^ £ \c T ajYa jS , c T Ts; S qq , J) (24) 
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TABLE I: Quantum numbers and masses for [cs][cs] states. We include in parentheses the mass shifts in MeV due to 
turning off the spin-spin interactions whence orbital angular momentum effectively increases the diquark-antidiquark 
distance. The dominant two-body decay channels are also indicated together with the relative angular momentum of the 
two produced particles and the partial width (F par t) of the exchange diagram dominated decays. As for the estimated 
masses, we put in parentheses the corrections in MeV due to turning off the spin-spin interactions whence orbital angular 
momentum effectively increases the diquark-antidiquark distance. By the notation we mean D or D* . 
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FIG. 1: [cs][cs] spectrum and decay thresholds of prominent channels. Dashed lines are for P— wave decays, solid lines 
for S— wave decays. As commented in [J, with the parameters at hand, we predict a mass value for the Y(4260) of 
m Y = 4330 ± 70 MeV. 



and similarly for the other operators. 

Results and discussion. The final results on the mass spectrum determination are summarized in Tab. |T] 
together with the spin, orbital quantum numbers, decay modes and widths, when calculable. We include 
in parentheses the mass shifts in MeV due to turning off the spin-spin interactions whence orbital angular 
momentum effectively increases the diquark-antidiquark distance. For a pictorial summary see also Fig. [T] 

Suppose that the tetraquark system could be described energetically by a double-well potential with the two 
light quarks lying in the two wells induced by the charm quarks, which can be considered as static color sources. 
The potential barrier separating the two wells prevents a quark in the diquark to bind with an antiquark in 
the antidiquark (and vice- versa). This process occurs anyways at the rate of the barrier penetration. We shall 
assume that this is the case for the quark passing anyhow through the barrier to bind with the antiquark in the 
other well, as represented in the following diagram 
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where X ss — <j>, u>, rj, rf. The charm quarks have no other choice than neutralize color in a charmonium meson. 
An alternative process is the formation of open charm mesons by rearranging of the strange quarks. We shall 
suppose that the latter two alternatives occur at almost the same rate (both of them pay the same energetic 
price of breaking the diquark bindings) . 

As for the decay widths, consider for example the decay into J/ip<fi °f the h particle with mass M = 
4277 MeV in Table H The S-matrix element is 

(JM^')0(e,<z)|W(p)) =g^ Vtl e uPp q a (25) 

where Q must have dimensions of the inverse of a mass to let the width have dimensions of energy. The decay 
at hand is a P— wave decay, therefore 

r(n 277 ^J/v^)4^(^ 3 ) (26) 

where p* is the decay momentum in the reaction and we assumed Q — g/My, the dimcnsionlcss g being written 
as g = A/ (\/2My). Assuming that the exchange diagrams in [c(u, d))[c(u, d)] states have the same amplitudes as 
the ones in [cs][cs] L 5j, i.e., assuming A = 2.6 GeV, we obtain that the decay modes reported in Table[U Similarly, 
S— wave modes, like the decays into J/^t? of the 1 H states in Table HI would be associated to a matrix element 
of the form (J / \p(r) , p')r](q)\Y (p, e)) = T where T must have dimensions of mass (one can set T = gMy). Matrix 
elements of the form (J/tl;(rj,p')r](q)\Y(p, e)) = N(p ■ T])(q ■ e) would instead give the D— wave contribution. As 
for P— wave decays like those of the 0~~ states we use the parameterization (J/i()(e, P ')ri(q)\0~~) = TL(p ■ e), 
where the dimensionless H is 7i = A/^/2. The standard 7777', ujtfi mixing schemes are used in the calculation. 

As another possible decay mechanism consider the quark pair creation allowing a diquark- antidiquark system 
to decay into a tetraquark and a standard meson as in [cs][cs] — » /o(980) J/ip. With quark pair production 
one could also have baryon-antibaryon final states. Charmed baryons with strangeness (S c , 0° c ) are anyway 
too heavy to be found in the decay products of the hadrons in Table fl] If we assume that the pair creation is 
regulated by a s (m c ) rather than by a s (AQCD), we might infer that pair creation is less probable. Along the 
same lines we do not expect to have significant contributions form 6-fermion interactions, as those induced by 
instantons (see Q). The reason is that the instanton Lagrangian is dumped by a factor exp(— 8n/g 2 ) which 
turns out to be rather small at the mass scale of the charm quark where a s {m c ) ~ 0.3. In this sense the 
instanton interactions are mainly related to infrared physics. Annihilation diagrams, also expected to be rather 
suppressed at the charm quark scale, could produce final states as DDtt. 

Experimental evidences. Invariant mass spectra of several of the final states of interest have already been 
explored by experimental searches. 

The most interesting match between predicted and observed states is in the J/ipu invariant mass spectrum, 
studied both by Belle [8j and by BaBar [9| in B — > J/tpuiK decays. A state with mass my = 3913 ± 4 MeV, 
according to the more accurate BaBar measurement, is observed to decay predominantly in this final state. 
This paper shows how the J PC = ++ state, which decays predominantly into J/ipui and can be produced in 
B decays in pair with kaons if L = 0, has a predicted mass of 3927 MeV. 

The J/tpr] invariant mass was studied by BaBar [l(| in B — * J/tprjK decays, and the resulting background 
subtracted distribution is reported in Fig. [2j Vertical lines refer to predicted mass value for the J PC = 
states. Even if states with different J PC quantum number may decay in J/tpn final state, selection rules forbid 
J PC = Q— or J PC = 1+- states. 

The V'/o(980) decay mode was instead studied when /o(980) decays into two pions, i.e. in the ^tttt final state, 
where ip can be either J/t/j or ip(2S). Exotic mesons are searched in Initial State Radiation (i.e. e + e~ — > V7 
processes) and can therefore only be J PC = The published spectra [HI, [13] show several structures at 

m = 4260, 4350, and 4660 MeV. Although only the latter one shows irir invariant masses clearly consistent with 
an /o(980) production, it is interesting to notice that the invariant masses predicted here for J PC = states 
are in the same mass range. 

Finally, possible exotic states decaying into J/ijxp and J/tpT)' have been searched in Refs. [T^] and [l4[ respec- 
tively but no significant signal has been observed even integrating over the mass spectrum. 
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FIG. 2: J/tpri events at BaBar. 



Conclusions. In this paper we have studied the consequences of allowing [cs][cs] diquark-antidiquark par- 
ticles with different J PC quantum numbers. We present their spectrum and main decay modes in Table [U 
The comparison with existing data shows some hints of match between observed and predicted particles but 
significantly larger data are needed for conclusive statements. 

After our paper appeared, the CDF collaboration presented a 3.8 a evidence of a narrow structure in J/ip4> 
at about 4143 MeV [15|. As it is clear from Fig. ® we do not have anything close in our spectrum. Other 
interpretations of the CDF structure can be found in [16j. In the same CDF paper, another peak in J/ip(f> is 
found with a significance < 3 a at about 4277 MeV. In this case we predict a h state at 4277 MeV decaying 
into J/ijxj> m -P— wave. 
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